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ABSTRACT 


The  principal  objective  of  this  thesis  project  was 
to  design  a  practical  remote  error  signal  generator  that 
could  be  incorporated  into  an  automatic  control  system 
for  a  farm  tractor.  The  practicality  of  the  system  was 
related  to  the  economic  situation  of  the  day.  Since  this 
thesis  was  prepared  at  a  time  when  prices  were  rising 
rapidly  due  to  inflationary  tendencies  in  the  economy,  and 
some  prices  were  dropping  fast  because  of  technological 
advances  in  the  electronics  industry,  the  estimated  price 
of  the  error  signal  generator  had  to  be  put  in  a  range  of 
twelve  to  eighteen  thousand  dollars. 

The  detection  system  used  a  method  of  triangulation 
to  locate  the  tractor.  Infrared  detectors  were  used  to 
indicate  when  the  tractor  was  in  the  detector  field.  The 
angular  position  of  the  detector;  was  related  through 
a  gear  train  to  the  position  of  the  shaft  of  a  digital 
shaft  position  encoder.  The  output  of  the  encoder  was 
a  digital  indication  of  the  angular  position  of  the  tractor 
relative  to  the  base  line  with  an  axis  of  rotation  passing 
vertically  through  the  detector.  A  digital  computer  was 
used  to  calculate  the  cotangents  of  the  angles  and  to 
perform  logical  operations  on  the  information  obtained. 

The  control  system  was  found  to  be  stable  for  a 
range  of  velocity  commonly  used  in  field  operations.  'Bith 
greater  demand  put  on  the  detection  system,  higher  veloci¬ 
ties  could  be  attained. 
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CHAPTER  I 


INTRODUCTION  AND  THE  PROBLEM 

A  current  review  of  automation  in  agricultural 
field  operations,  shows  that  various  forms  of  mechani¬ 
zation  are  being  developed  to  compensate  for  the  short¬ 
age  of  labor.  The  emphasis  on  the  problem  will  un- 
doubtably  lead  to  the  development  of  automatic  control 
systems  because  of  the  significant  decrease  in  human 
involvement  possible.  This  statement  is  not  implying 
that  automatic  control  systems  do  not  already  exist  in 
agricultural  operations  but  that  the  sphere  of  applica¬ 
tion  is  limited  to  specific  minor  operations.  Major 
areas  exist  in  which  measurement  and  control  systems 
can  be  expected  to  play  an  increasing  role  in  the  future. 

Manufacturers  of  instrumental ion  and  control  equip¬ 
ment  are  ax\Tare  of  the  vast  market  for  their  services  in 
agricultural  production  and  are  anxious  to  exploit  this 
situation.  However  there  is  no  immediate  revolution  to 
be  expected  as  a  result  of  this  interest.  Nhen  problems 
of  the  environment  and  of  the  agricultural  production  have 
been  overcome,  there  still  remain  considerations  oi  price 
and  likely  demand  which  often  discourage  the  development 
of  commercial  equipment.  Nevertheless, 


the  continuance 


of  present  trends  in  agriculture  and  in  the  electronics 
industry  must,  in  time,  favor  further  developments.  The 
farmer's  growing  expenditure  on  mobile  and  static  equipment 
is  an  obvious  trend  of  this  kind.  Ensuing  economy  in  man¬ 
power  and  effort,  improvement  in  production  equality  and 
quantity,  and  further  use  of  favorable  weather  and  soil 
conditions  .  justify  considerable  outlay.  In  electronics 
an  encouraging  sign  is  the  growing  availability  of  the 
packaged  solid  state  units  and  integrated  circuits  for 
measurement  and  control  systems.  3y  employing  these  in 
agricultural  applications,  maximum  price  advantage  can  be 
attained  in  addition  to  simplified  construction  and  re¬ 
placement  . 

1.2  LOCATION  AND  CONTROL  OF  FIELD  TRACTORS 

The  general  purpose  tractor  is  the  main  source  of 
power  in  North  American  agricultural  field  crops.  Cur¬ 
rently  these  are  not  designed  to  give  the  driver  a  rea¬ 
sonable  degree  of  comfort  or  safety.  For  this  and  the  rea¬ 
sons  given  in  Section  1.1,  engineers  in  various  parts  of 
the  world  have  been  exploring  ways  to  dispense  with  the 
driver  for  some  types  of  field  operations.  Several  al¬ 
ternate  approaches  have  been  taken  in  attempts  to  limit 
the  level  of  human  intervention  in  the  control  of  the 
tractor. 

The  most  widely  explored  method  is  the  use  of  radio 
1  ? 

control  (  the  tractor  responding  to  commands  which 

(2) 


operate  tuned  relays  linked  by  electrical  or  electrohy- 
draulic  devices  to  the  controls.  Immediate  applications 
are  for  operating  the  tractor  in  areas  where  some  form 
of  hazard  exists,  and  control  of  more  than  one  tractor  by 
a  single  operator,  i,e.,  tandem  cultivation  with  the  oper 
ator  on  one  tractor  and  the  slave  tractor  following  close 
ly3.  For  more  remote  control,  closed  circuit  television 
is  necessary. 

The  radio  control  system  may  in  these  ways  increase 
output  or  make  the  operator’s  life  more  pleasant,  but  it 
does  not  offer  the  best  approach  to  a  more  ambitious  pro¬ 
ject  -a  fully  automatic  tractor  guidance  system.  The  at¬ 
traction  of  such  a  sytem  is  that  it  may  conserve  labor 
by  freeing  the  driver  -at  least  in  part-  from  the  more 
routine  operations  and  enable  him  to  carry  on  with  work 
requiring  a  greater  degree  of  human  intervention.  The  re 
gular  shape  and  topography  of  the  fields  in  Western 
Canada,  permit  a  number  of  possible  guidance  systems. 

A  direct  line  of  sight  system  employing  the  laser 
beam  is  described  by  Williamson The  drift  limit 
with  the  laser  detector  is  5/8  inch  in  2  miles.  A 
gas  laser  is  mounted  parallel  to  the  telescope  of  a 
surveyor’s  transit.  The  laser’s  narrow  continuous  beam 
is  aligned  by  the  transit  to  fall  on  the  center  oi  a 
rectangular  array  of  photocells  mounted  at  front  and 
back  of  the  machine.  The  photocells  are  tied  in  with 

(3) 


an  identical  array  of  lamps  on  a  panel  in  front  of  the 
driver.  Thus  the  illuminated  photocells  light  equivalent 
lamps  on  the  panel,  giving  the  operator  an  identication 
of  how  closely  he  is  tracking  the  laser  beam. 

At  the  NIAE^,  a  system  employing  ground  markers,  has 
been  investigated.  A  working  model  has  been  constructed 
on  a  tractor  with  an  infinitly  variable  hydrostatic  trans¬ 
mission.  A  similar  system  using  furrow  following  links 
is  described  in  Farm  Mechanization^.  Another  version  of 
ground  marker  is  the  one  using  a  leader  cable.  This  was 

developed  by  Morgan'7  at  Reading  University,  Finn-Kelcy8 
8  s 

and  Owen  ^  recently  presented  a  paper  on  an  automatic  sys¬ 
tem  of  programmed  tractor  control  using  the  same  technique 
as  Morgan. 

Another  dead --reckoning  control  system  which  was  fit¬ 
ted  to  a  tractor  was  described  by  Gilmour^,  Richey-10 
developed  a  system  that  operated  on  the  principle  of 
plant  feelers.  Although  this  system  did  not  eliminate 
human  intervention,  it  definitely  decreased  the  degree  of 

human  concentration  necessary. 

11 

MacHardy  described  an  automatic  guidance  system 
for  a  farm  tractor.  The  system  would  keep  the  tractor  on 
a  course  parallel  to  a  given  reference  line,  with  only  three 
level  (-,+,0)  position  information.  He  proposed  a  method 
for  determining  the  position  of  the  tractor,  but  did 


. 


not  research  this  aspect  of  the  project.  Shukla1^'.  at 
Missouri  State  University,  is  working  on  a  control  sys¬ 
tem  that  is  based  on  somewhat  the  same  principle,  but 
to  date  there  is  no  record  of  a  publication. 

1.3  OBJECTIVES  OF  THIS  THESIS 

The  objectives  of  this  thesis  are  as  follow: 

a)  To  describe  a  practical  system  for  determining 
the  position  of  the  tractor  in  an  automatic  con¬ 
trol  system  for  a  farm  tractor. 

b)  To  evaluate  the  precision  obtainable  with  the 
system. 

c)  To  estimate  the  cost  of  the  hardware  involved 
in  the  system. 

1.4  OUTLINE  OF  THE  PROJECT 

A  synopsis  of  the  automatic  guidance  system  to 
which  the  error  detection  system  Is  to  be  applied  will 
follow  in  the  next  chapter.  The  characteristics  re¬ 
quired  in  the  error  generator  will  then  follow.  The 
The  equipment  will  be  described  generally  and  mathe¬ 
matically  in  ensuing  sections.  The  last  sections  will 
be  an  evaluation  of  the  response  of  the  complete  system, 
to  a  set  of  initial  conditions. 
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CHAPTER  II 


SYNOPSIS  OP  THE  GUIDANCE  SYSTEM 

A  tractor  guidance  principle  that  served  as  a  starting 
point  for  the  research  leading  to  this  thesis,  is  described 
in  a  paper  published  by  MacHardy11. 

MacHardy  simulated  the  tractor  control  system  on  an 
analog  computer,  and  built  a  working  tractor  model. 

In  the  paper  mentioned  above,  MacHardy  also  proposed 
a  method  for  locating  the  tractor  by  means  of  infrared 
detectors.  The  system-  for  locating  the  tractor  was  not 
investigated  however.  The  extension  of  the  project  to 
include  the  system  for  locating  the  tractor  is  the  prime 
object  of  this  thesis, 

2.2  REVIEW  OF  THE  GUIDANCE  SYSTEM  THEORY 

A  review  of  the  control  system  referred  to  (11,13) 
follows . 

MacHardy  proposed  to  guide  a  tractor  along  the  path 
rr’ ,  Pig.  2.1,  which  is  parallel  to  the  reference  line  AB. 


-I _ _  2?  1 


FIGURE  2.1  LOCATING  THE  POINT  P  BY  TRIANGULATION 
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_ 
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The  distance  between  the  observation  points  A 
and  B  is  known.  ?  is  the  actual  position  of  the  tractor; 
therefore  it  is  a  distance  h  above  the  baseline  AB. 

Since  the  sum  of  cot^c  and  cot£  is  constant 
a/h+b/h  =  A3/h 

for  all  points  on  a  path  parallel  to  A3, 


A3 


=  h 


2.2.1 


cot^  +  cot  £ 

Should  the  distance  h  be  different  from  any  desired 
tractor  location  ,  an  error  can  be  measured  as 

e  =  AB/(cot~  +  cot?  )J  -  2.2.2 

NacHardy  proposed  that  the  error  e  should  have 
three  levels,  (+,-,0)  depending  on  the  location  of  the 
tractor. 


m 

.e  n  0 

0- 

8u 

» 

FIGURE  2  o 2  ERROR  SIGNAL  OUTPUT  AS  A  FUNCTION 

OF  THE  DISPLACEMENT  y 

Figure  2.2  illustrates  the  error  signal  as  a 
nonlinear  function  of  the  lateral  tractor  displacement 
e  where  eu  and  e^  represent  the  upper  and  lower  limits 
within  which  it  is  proposed  to  control  the  tractor. 

This  signal  determines  the  tractor  front  wheel  angular 
position  as  turned  to  the  left,  straight  ahead  and 
turned  to  the  right.  The  tractor  response  to  a  turning 

(7) 
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signal  may  be  described  as  follows: 


Let  Tm  be  the  anSular  position  of  the  front  wheels 
in  radians  per  foot,  and  let  Ta  be  the  angular  position 
of  the  tractor  with  respect  to  the  baseline  in  radians. 

Th  en 

t 

Ta  =  /  T«dt  2.2.3 

where  v  is  the  forward  velocity  of  the  tractor  in  feet 

per  second.  Let  y  be  the  lateral  distance  from  the  refer¬ 
ence  line.  Then 


y  -  y/  sinTadt 

For  Ta40.20,  sinTa  -  Ta < . Ols in( Ta) ; 
hence  for  small  angles,  let  sin(Ta) 
Th  en 

tt 


y  -  vyf  Tmdtdt 


00 


2.2.4 


■a' 


2.2.5 


It  was  proposed  to  damp  the  steering  system  by  using  a 
gyrocompass  in  a  minor  feedback  loop.  The  gyrocompass 
would  return  some  part  of  the  velocity  with  which  the 
lateral  position  of  the  tractor  was  changing.  If  the  damp¬ 
ing  is  set  to  b,  then  the  equation  developed  by  MacHardy, 
2.2.5,  describing  the  steering  control  system,  becomes 


y  =  T  e“Dtdt 
0 

2.3  GENERAL  ANALYSIS 
If  the  open  loop 
system,  Eq.  2.2.6,  is 

,  ,  c(s) 

g(s)  =  -Hr(sj  = 


2.2.6 

OF  THE  INTEGRO-DIFFERENTIAL  EQUATION 
transfer  function  of  the  guidance 
expressed  in  the  Laplace  domain;  it 
Tmv^ 

~(s+  rj 
(8) 


becomes 


2.3,1 


■  ; 


in  which  Tm  is  a  coefficient  that  can  have  only  one  of 
three  discrete  values  at  a  time.  The  error  signal 
E(S)  =  R(S)  -C (S )  2.3  2 

Eliminating  E(S)  from  Eq.  2.3.1  and  letting  K  = 
results  in 


(S(S+1))C(S)  +  KC(S)  =  KR(S) 

Taking  the  inverse  Laplace  transform. 


j~-  C  ^  c  (  ^ )  .{-  Kc  ( t )  ~  Kr  ( t ) 

dt2  dt 


If  r(t)  is  zero,  Eq. 
2 

d  c(t)  +  do ( t )  + 
dt2  dt 


2. 3.4  becomes 
Kc  ( t )  =  0 


2.3.3 


Eq  2.3.3  becomes 
2.3.4 


2.3.4 


Substitute  x=c(t)  and  y=  £  in  Eq.  2.3.5.  Then  the 


normalized  equation  becomes 
y  =  ~y  -Kx 

Divide  Eq.  2,3.6  by  v  =  x;  then 


2.3.6 


The  last  equation  presents  an  analytical  expression  for 
the  slope  ox  a  trajectory  In  the  (x,x)  plane.  For  an 
input  r(t)  =  0,  there  exists  one  trajectory  in  the  phase 
plane  for  each  set  of  Initial  conditions  x(0),  and  x(0). 

Absence  of  any  forcing  function  in  Eq.  2.3.?  implies 
that  all  sytem  excitation  exists  in  the  form  of  an  init¬ 
ial  lateral  dislocation  and  angular  position  of  the  trac¬ 
tor.  In  the  phase  plane  the  trajectory  starts  at  the  point 
described  by  the  initial  conditions.  Using  isoclines  (loci 
In  the  phase  plane  that  have  constant  slope),  the  plot 
a  trajectory  in  the  phase  plane.  Fig.  2.3  is  obtained. 
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Table  I  contains  information  on  the  slope  of  the  trajec¬ 
tory  in  the  region  of  the  isoclines  in  the  phase  plane. 
These  regions  of  constant  slope,  facilitate  plotting  the 
complete  response  of  the  system  for  a  set  of  initial 
conditions.  The  initial  conditions  of  the  trajectory. 
Fig.  2.3*  &re  (4.0,  0);  in  which  the  4  represents  a  lat¬ 
eral  displacement  of  the  tractor  of  4  feet  and  the  0 
indicates  that  the  lateral  position  of  the  tractor  is 
not  changing. 


TABLE  I 


TA3LE  OF  VALUES  SATISFYING  dy/dx  =  -1  -Kx/y 


dy/dx 


y  =  f(x) 


-6 

“5 

-4 

0 

1 

2 

3 

5 

6 


1 .  Ox 

1  o  3x 

408x 

-5*9  x 
-2 ,8x 

-1  o8x 
~1 . 4x 
~0o9x 

-0.7x 


y 


X 
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2.3,  demonstrates 


The  phase  plane  trajectory,  Fig. 
graphically  that  the  system  described  by  MacHardy  was 
stable  for  some  initial  conditions. 


i 
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CHAPTER  III 


DETECTION  SYSTEM  COMPONENTS 


Nhereas  the  choice  of  components  for  the  guidance 
system  developed  as  alternatives  were  considered,  (and 
often  rejected),  it  is  felt  that  there  should  be  a  brief 
description  of  the  system’s  components  as  finally  de¬ 
cided  upon,  before  discussing  the  reasoning  that  lead 
to  the  particular  choices. 

Each  infrared  detector  is  scanned  horizontally 
through  a  small  arc  about  the  target  position.  The  an¬ 
gular  position  of  each  detector  is  measured  by  a  shaft 
position  encoder  with  a  digital  output.  The  encoder 
output  is  connected  to  interfacing  equipment  and  a 
PDP-8  digital  computer.  The  computer  is  programmed  to 
compute  the  cotangents  of  the  angles  locating  the  tar¬ 
get  from  both  detectors,  then  to  sum  these  cotangents 
and  compute the  tractor  lateral  location  as  indicated 
in  Section  2  03*  The  actual  lateral  position  is  com¬ 
pared  with  the  programmed  position,  and  if  these  are 
unequal.,  an  error  signal  is  radioed  to  the  tractor 
steering  mechanisim  Figure  3»1  shows  a  simplified 

diagram  of  this  system  and  its  componentso 

-From  second 
'detector  and 

encoder 


Infrared - ^-Detector 

radTati\m_^ _ v 

fab* 


r 


purer' 

Radio 

Compute 

r 

X-mit 

Encoder 


Steering  signal 


to  tracrbor 
FIGURE  3.1  ERROR  DETECTION  SYSTEM 
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3.2  GENERAL  DESCRIPTION  0?  THE  DETECTOR 


Infrared  detectors  are  sensing  devices  that  de¬ 
tect  objects  whose  surface  temperatures  are  effectively 
different  from  that  of  their  surroundings.  The  first 
concentrated  interest  in  their  development  was  evi¬ 
dent  during  the  first  world  was  when  the  German  scien¬ 
tists  began  exploring  their  use  in  military  applica¬ 
tions.  During  the  second  world  war,  the  interest  again 
increased  and  devices  for  following  equipment  movements 
at  night  were  developed  and  improved. 

Today  the  military  applications  of  detectors  still 
outnumber  the  industrial  and  medical  applications. 
Missle  seekers  and  homers  (used  in  Sidewinder  air  to 
air  missies)  J  use  infrared  detectors  to  find  and  lock 
on  hot  areas  of  the  intended  target.  Other  military 
applications  are  fire  control,  bomber  defence,  and  the 
ballistic  missle  detection,,  The  two  main  reasons  for 
their  popularity  in  military  applications  are  (l)  the 
passive  nature  of  their  detection  system  and  (2)  the 
increased  angular  refinement.  Unlike  radar,  they  usu¬ 
ally  do  not  emit  a  signal.  This  renders  them  undetec¬ 
table.  Active  systems-^  do  exist;  however  they  are  not 
widely  used. 


Peaceful  applications  are  widely  diversified. 

One  detector  developed  by  the  NRG  is  used  to  detect  hot 
spots  in  high  voltage  transmission  lines.  The  forestry 
departments  use  detectors  to  locate  accurately,  the 
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fires  shrouded  by  heavy  smoke,  thus  increasing  the 
effectiveness  of  the  fire  fighting  equipment.  The 
detector  that  was  developed  by  the  NRG,  is  the  one 
after  which  the  detector  that  was  used  in  the  ex¬ 
perimental  work  leading  to  this  thesis, was  designed. 

There  are  many  more  applications  and  a  few  other 
types  of  detectors  in  use;  however  mention  has  been  re 
stricted  to  those  applications  that  are  closely  re¬ 
lated  to  the  one  used  in  this  thesis. 

The  active  components  of  the  Bolometer  type  in¬ 
frared  detector  consist  of  thin  resistive  elements  ex¬ 
hibiting  a  large  change  of  resistance  with  temperature 
connected  in  a  bridge  circuit  as  in  Pig.  3*2° 


FIGURE  3  o  2  BOLOMETER  BRIDGE  CIRCUIT  ARRANGEMENT 
The  sensing  element  in  the  detector,  is  a  semi¬ 
conductor  film  of  a  mixture  of  manganese,  nickel  and 
cobalt.  These  thermistors  have  Tar gel  negative  tem¬ 
perature  coefficients.  A  rapid  response  is  obtained 
by  attaching  the  "flakes"  to  good  heat  conducting 
thermal  sinks. 

The  bolometer  bridge  circuit,  Pig,  3*2,  comprises 
two  identical  thermister  elements  mounted  on  the  same 
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base  to  effect  compensation  for  changes  in  ambient  tem¬ 
perature.  One  element,  the  active  reciever,  is  exposed 
to  infrared  rays  while  the  other  is  shielded. 

3.3  INPUT -OUTPUT  RELATIONSHIP  OF  THE  DETECTOR 

Fig.  3.2  (a)  shows  the  functional  units  of  the 
electronic  system. 


,A 


(a) 


(b) 


FIGURE.  3. 3  (a)  FUNCTIONAL  UNITS  ARRANGED  IN  A  BRIDGE 
CIRCUIT;  (b)  THEVENIN  EQUIVALENT  CIRCUIT. 

When  the  voltage  source  is  set  to  zero,  the  equiva¬ 
lent  impedance  seen  at  terminals  AB  consists  of  the 
parallel  combination  of  Z£  and  z^ ,  in  series  with  the 
parallel  combination  of  z^  and ' .  Thus  in  equation  form 


Z  «  =  z 1 z 4  +  Z2Z3 

zi  '**  z4  z2  +  z3 


3.3.1 


On  open  circuit  the  source  results  in  currents 
1^  and  I 2  as  shown  in  Fig.  3*3* (a). 

V 


Ii  = 


VS 

z^+  z  4, 


and  I 2  = 


cr 


3.3*2 


z2  1*  z3 

Assuming  the  potential  of  A  greater  than  B, 


V 


t  — 


VAB  =  Ixz 4  -  I2z3 
= 


Zi  +  Z ,+ 


Zo  +  Z 


3.3.3 

3.3.^ 


3 
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V* 


3.3.2 


=  V, 


Z2Z4  -  z  1 z  3 


*  (z!  +  Z4)(Z2  +  Z3) 

The  Thevenin  equivalent  voltage  Vs 


the  difference  Z2Z4  “  zlz3 
voltage  V*  =  0 


.  When 


is  proportional 
z-^z^  =  z2z4»  the 


to 


3.4  TIMS  RESPONSE  OF  THE  DETECTOR 


The  level  of  response  that  the  detector  reaches  is 
dependent  upon  two  factors:  (l)  the  position  of  the  target 
in  the  detector  field  and  (2)  the  duration  that  the  target 
remains  in  the  field  of  the  detector.  The  response  of  the 
detector  as  a  function  of  the  target’s  position  in  the 
field  of  the  detector  is  illustrated  in  Fig.  3*4. 


FIGURE  3.4  RESPONSE  OF  THE  DETECTOR  AS 
OF  THE  TARGET  POSITION  IN  THE  DETECTOR 


A  FUNCTION 
FIELD 


The  response  curve  shown  in  Fig.  3.4  is  a  reproduc¬ 
tion  from  a  figure  prepared  by  personnel  at  the  raanu- 
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facturing  laboratory.  The  curve  can  be  approximated 
by  Eq.  3.4.1  in  which  the  magnitude  is  y;  c  is  a  con¬ 
stant  due  to  imbalance,  and  the  zero  point  on  the  axis 
is  at  the  center  of  the  response  curve.  The  width  of 
the  detector  field  is  denoted  by  r  and  the  factor  K  is 

(16) 


Eq.  3.4,1,  is  l/r/2  for  0<irr/2  and  0  for  ©>-rrr/2,  with 
expressed  in  radians. 


r/2 

y  =  c  +  2n/  -sinKGdG 
-r/2 


3.4.1 


r 


n  is  a  factor  that  is  itself  a  function  of  time.  The 

factor  n  is  very  near  linear  in  its  magnitude  with  respect 

to  time  where  the  time  does  not  exceed  2T~  where  T  is 

c  c 

the  time  required  for  the  detector  to  attain  a  level  of 
response  that  is  50/  of  its  final  response  to  a  partic- 
ular  disturbing  temperature  change. 

b 

n  =/n(t)dt  for  b<2Tc  3A.2 

n  =  nk  for  b?2Tc  3.4.3 

There  need  be  no  concern  with  n  in  Eq.  3. 4.2  if  the  width 
of  the  field  w  and  the  frequency  u>  are  such  that  the 
detector  will  be  on  target  for  time  t>  2T0.  Thus  the 
position  of  the  detector  may  be  described  by 

0  =  w/2  x  (sin^t)  3.4.4 

To  obtain  the  50/  level  of  response  from  the  detector 
at  a  particular  angular ‘ position  with  respect  to  the  tar¬ 
get  source,  w  and  a;  may  be  varied.  Thus  for  example,  if 
it  is  desira ole  to  have  the  50/  level  of  response  occur 
at  the  moment  the  tractor  is  on  the  center  of  the  deted- 
toi  field,  the  angular  velocity  of  the  detector  scan  is 
made  such  that  the  response  would  reach  50/  at  that  time. 

If  the  frequency  of  scan  uj-  is  increased,  there  must  be 
a  corresponding  decrease  in  the  angular  width  of  the  field 
scanned  in  order  to  maintain  the  50/  response  at  the  same 
relative  position  of  the  tractor  in  the  detector  field. 

(1?) 


3.5  FIELD  OF  VIEW  AND  SENSITIVITY  OF  THE  DETECTOR 

The  field  of  view  is  defined  as  the  angular  width 
by  height  of  the  region  of  sensitivity.  The  sensitivity 
is  defined  as  the  ability  of  the  detector  to  distinguish 
the  target  radiation  from  the  infrared  rays  emmitted 
by  the  background. 

Since  a  tractor  could  only  be  located  accurately 
if  the  detector  could  distinguish  between  the  exhaust 
pipe  and  the  cloud  of  exhaust  gases,  it  was  felt  nec¬ 
essary  to  test  the  detector* s  ability  to  distinguish 
between  these  sources  of  heat.  Trials  with  the  detec¬ 
tor  were  conducted  by  the  author  in  1966.  The  purpose 
of  the  trials  was  to  determine  the  practicality  of  using 
this  mode  of*  detection.  The  direction  of  focus  for 
the  tests  was  north.  The  reason  being  that  the  highest 
degree  of  interference  by  the  background  is  antici¬ 
pated  to  arise  from  the  south  side  of  objects  in  the 
northern  hemisphere.  Another  purely  mechanical  reason 
for  the  choice  of  direction,  was  to  prevent  accidental 
exposure  of  the  lens  to  direct  rays  of  the  sun. 

In  the  trials,  a  simulated  tractor  exhaust  muffler 
was  used.  The  simulated  muffler  was  a  rusty  heating 
jacket  and  fire  ring  from  a  side  arm  heater.  A  variable 
speed  motor  with  a  centrifugal  fan,  forced  air  through 
the  heater.  Propane  was  then  adjusted  until  an  exhaust 
gas  temperature  comparable  to  that  of  a  tractor  under 
loaded  conditions  was  attained.  This  meant  that  there 
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was  an  excess  of  air  of  approximately  200%  of  that 
required  for  complete  combustion.  Data  to  show  the 
effect  of  the  exhaust  gases  on  the  detector  response, 
are  contained  in  Table  II, 


TABLE  II 


OUTPUT  OF  THE  DETECTOR  AMPLIFIER  WITH 
VAPORS  AND  HEATER  EXPOSED,  WITH 
HEATER  EXPOSED  AMD  WITH  JUST 
THE  VAPORS  EXPOSED 


Distance 

yards 

Output  i 
micro-’ 
amps 

attenuator 

setting 

Direction 
of  view 

Field 

width 

Background 

220 

2.5 

.5 

north 

0  °  1 2  ’ 

trees 

220 

12.0 

1.0 

tt 

0°12» 

tt 

Heate 

r  and  vapors  exposed 

220 

2.5 

.3 

north 

0°12 ' 

n 

220 

20,5 

1,0 

n 

0°12  1 

1! 

Heater  expos 

ed  only.  Vapors  shielded  by- 

a 

4x8  piece 

of  ■§■  in. 

plywood . 

220 

0 

>3 

ii 

0°X2 ' 

It 

220 

•  5 

1.0 

0°12» 

tt 

Vapors  plus  i  in.  of  the  heater  were  exposed. 

In  a  sweep  of  3§0o  the  variation  of  the  meter  was  . 2j;.amp 


Before  the  above  results  were  obtained,  it  was 
necessary  to  prevent  fluctuation  of  the  temperature 
of  the  detector  barrel,  especially  in  the  vicinity  of 
the  bolometers.  The  temperature  fluctuation  was  mini- 
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raized  by  shielding  the  barrel  from  the  wind  with  a 
layer  of  insulating  material. 

The  data  in  Table  II  indicate  that  there  is  an 
insignif icant  contribution  by  the  vapors  to  the  total 
infrared  radiation  received  by  the  detector.  When 
the  vapors  were  shielded,  the  plywood  was  heated 
sufficiently  to  cause  the  total  width  of  the  sensitive 
region  to  be  increased. 


TABLE  III 


SENSITIVITY  AT  A  DISTANCE  0?  ONE  QUARTER  MILS 


Distance 

yards 

Output 

micro- 

amps 

Attenuator 

setting 

Direction 
of  view 

Field 

width 

Background 

440 

.5 

•  3 

north 

j 

CO 

O 

1  0 

trees  and 
cattle 

440 

4o  0 

1.0 

m 

CO 

O 

O 

Temp  60° 

440 

2.5 

1.0 

11 

CO 

0 

0 

Temp  ?0°  = 

Width  of  target  0°2* 

Heater  and  vapors  were  exposed  for  these  and  the 
remaining  trials. 

Less  than  0.5  microampere  variation  in  a  sweep 
of  360°. 


Comparing  the  data  of  Table  III  with  that  of  Table  II 


confirms  the  law  that  the  intensity  of  infrared  rays 
decreases  as  the  square  of  the  distance.  However  there 
is  still  substantial  output  from  the  detector  amplifier 
at  a  distance  of  one  quarter  mile. 

In  Table  IV  it  is  seen  that  the  ability  of  the  de¬ 
tector  to  distinguish  the  target  under  windy  conditions 
was  no  different  from  that  on  still  days.  Stability  of 


(20) 


the  detector  system  was  substantially  improved  by  preven- 


TABLE  IV 

SENSITIVITY  UNDER  WINDY  CONDITIONS 


.stance 

Output 

Attenuator 

yards 

Micro- 

Setting 

amps 

440 

Variable 

i.O  The  detector  amplifier 

440 

205 

response  Is  unstable 
1.0  Barrel  now  protected 

from  the  wind  by  styrofoam 
insulating  material. 


ting  fluctuation  in  temperature  of  the  detector  barrel 
especially  in  the  vicinity  of  the  bolometers. 

TABLE  V 


DETECTOR  AMPLIFIER  OUTPUT  WITH  CHANGING  DIRECTION 


Distance 

yards 

Output 
micro - 
amps 

Attenuator 

setting 

Direction 
of  view 

Position 
of  sun 

Tem¬ 

pera¬ 

ture 

440 

2.5 

1.0 

north 

south 

6  0° 

440 

2.5 

1,0 

east 

s6o°s 

50° 

440 

4.0 

1.0 

N45°E 

S55°3 

52° 

The  background  for  the  first  two  trials  was  trees; 

for  the  last  trial  above,  the  background  was  sky. 

Note  the  substantially  higher  definition  of  the  target 
with  a  background  of  clear  sky.  The  direction  of  focus 
relative  to  the  sun  has  apparently  little  or  no  significant 
effect  on  the  total  infrared  radiation  of  the  background. 

Table  VI  indicates  that  the  sensitivity  of  the  detec¬ 
tor  improved  substantially  in  darkness.  A  surprising 
observation  is  that  during  a  period  of  moderate  rain,  the 


(21) 


TABLE  VI 

SENSITIVITY  VARIATION  WITH  DARKNESS 


Time  Output 
micro- 
amps 

Attenuat or 
setting 

Direction 
of  view 

Background  Weather 

4PM  3.0 

1.0 

north 

barley 

rain  70° 

8PM  5.0 

1.0 

ti 

it 

"  64° 

12:30AM  6.0 

1.0 

n 

ii 

no  rain 

65° 

7:45AM  5o0 

1.0 

it 

it 

"  »  68° 

2:00PM  4.0 

1.0 

it 

it 

11  it  75^ 

No  varia 

tion  detectable  for  360° 

sweep  at 

night , 

sensitivity  was  not  noticably  decreased. 

3.6  TIME  CONSTANT  OP  THE  DETECTOR 

To  check  the  time  constant  of  the  detector,  ap¬ 
paratus  was  arranged  as  in  the  block  diagram,  Fig.  3.5* 


FIGURE  3,5  BLOCK  DIAGRAM  OF  THE  APPARATUS  USED 
TO  CHECK  THE  TIME  CONSTANT  OF  THE  DETECTOR 

The  camera  lens  was  removed  and  the  shutter  was  held 

open  while  the  lamp  was  adjusted  until  a  steady  output 

was  reached  on  the  plotter.  The  shutter  was  then  closed 

to  plot  the  zero  position.  The  time  base  function  was 

produced  by  integrating  a  constant  with  respect  to  time 

and  applying  this  signal  to  the  X-axis  of  an  X-Y  plot¬ 


ter. 
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FIGURE  3.6  TIME  RESPONSE  OF  THE  DETECTOR 
The  response  curve  of  the  detector,  Fig.  3.6,  confirms 
that  the  time  constant  of  the  detector  is  very  close  to 
the  manufacturers  specification  of  6.5  milliseconds.  The 
6.5  millisecond  level  is  about  il  millivolts  lower  than  the 
the  50 %  level  of  the  steady-state  response however  the 
discrepancy  can  be  charged  partly  to  the  plotter. 
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3.7  SPECIFICATIONS  CONTINUED:  FIELD  OF  VIEW 

The  manufacturer’s  specifications  lists  the  field, 
of  view  of  the  detector  as  0.1°  wide  x  1.5°  high. 

Trials  conducted  with  the  instrument  by  the  author, 
(Section  3*5)  confirmed  this  specification.  The  angular 
resolution  of  the  detector  is  affected  by  the  width 
of  the  field  of  view. 

For  example,  consider  a  tractor  as  in  Fig.  3»7 
with  the  detectors  located  at  distances  A  and  3  from 
the  tractor,  the  widths  of  the  detector  fields  are 


Sa  and  S-q  respectively. 

si  =  ^i®  ~  VTi  (  i  -  a  or  b) 


3.7.1 


in  which  9  is  in  radians  and  w  is  the  width  of  the 
field  of  view  at  a  distance  from  the  detector  in 
feet . 

Sa  =  1320  x  .001164  =  1.536  ft o 
sb  =  1865  X  .001164  =  2.173  ft. 

The  area  in  which  the  tractor  is  located  is  a 
parallelogram  with  sides  of  1.536  ft.  and  2.173  ft. 
as  shown  in  Fig.  3*7.  The  parallelogram  has  a  width 
of  4.50  feet  perpendicular  to  the  intended  tractor 
path  rr’ .  This  exceeds  the  desired  limit  of  plus  or 
minus  one  foot,  thus  further  refinements  in  the  system 
are  necessary  in  order  to  achieve  the  desired  limit. 
This  is  referred  to  again  at  a  later  point. 
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FIGURE  3.7  TRACTOR  EXHAUST  LOCATION  LIMITS 


3.8  DETECTOR  SCANNING  VELOCITY 

The  time  constant  of  the  detector  sets  an  upper 
limit  on  the  angular  velocity  of  the  scanning  system. 
In  the  absence  of  any  auxiliary  scanning  devices,  the 
time  constant  of  the  detector  will  not  permit  a  contin 
uous  sweep  of  the  field,  thus  leading  to  the  oscil¬ 
lating  drive  system  that  is  represented  in  the  block 
diagram  of  Fig.  3*8* 


FIGURE  3.8  BLOCK  DIAGRAM  OF  THE  DETECTOR  DRIVE 

SYSTEM 

In  the  aljove  system,  the  feedback  for  reversing 
the  motor  is  provided  by  the  detector.  The  detector 
feedback  is  input  to  a  Schmitt  trigger  and  the  output 
from  the  Schmitt  trigger  controls  the  power  in  the 
drive  circuit.  This  system  ensures  that  the  detector 
will  follow  the  tractor,  while  sweeping  the  region 
in  which  the  tractor  is  located  with  an  oscillating 
motion . 
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Under  normal  operating  conditions  the  velocity  of 
scan  would  be  limited  to  0.22  radians/second.  At  this 
velocity  and  considering  the  time  delay  due  to  the  de¬ 
tector  time  constant,  the  detector  will  produce  its 
pulse  when  the  tractor  is  at  the  center  of  the  detector 
field.  Thus  a  higher  degree  of  precision  is  achieved. 


3.9  LATERAL  FREEDOM  DUS  TO  NON SYNCHRONIZATION  OF  THE 
DETECTORS 

The  two  detectors  in  the  system  result  in  an 
added  problem. 

If  the  tractor  has  a  velocity  of  v  feet/second 
and  the  detectors  are  out  of  synchronization  by  t 
seconds,  the  sum  of  the  cotangents  will  be 

a/h  +  (A3  -  a  ±  vt0)/h  or  3*9*1 


cot  (a  -  )  +  cot  (in  +  tQAt)  =  AB  +  vtQ  3*9*2 


h 


This  can  also  be  written  as 


h  =1 


A3  +  vt, 


AB  -vt 


o 


cot ( a;jL)*i-cot  ( b At 0 )  cot  ( ajL)+cot ( bj_-Afc0  j 


3*9*3 


For  example,  suppose  that  the  tractor  is  being  operated 
with  velocity  8  feet/second  at  point  P,  660  feet  from 
the  base  line  AB  in  Fig.  3*9 »  when  the  first  detector 
(a)  registers  an  output  signal.  T0  seconds  later  de¬ 
tector  (b)  registers  its  signal. 


(2?) 


* 


. 


Then 


cot(a)  -  660/660  =  1 
cot(b)  =  (660  +  8)/66o 


Table  VIII  gives  the  lateral  freedom  for  various 


time  lags  tQ. 


TA  p  ttt 

-i.  2*  ID  a— <— 1 


v  . 


:ii 


LATERAL  FREEDOM  OF  Tr 
S YN CHRON IZAT I ON  0 F 


S  TRACTOR  DUE  TO  NON- 
THE  DETECTORS. 


Time  Sum  of 
period  Cotangents 
sec . 


AB/c  o  t  ( a )  -r  c  ot  ( b ) 


1 

2 


Lateral 

freedom 

feet 


(hi-tc  -hi+ 


o 


X 

) 


.0 

2.0121212 

1.9878787 

656.024 

664*  024 

4.024 

.75 

2.0090909 

1.9909090 

657.013 

663.013 

3.013 

.50 

2.0020202 

1.9939393 

658.006 

662 . 006 

2.006 

.25 

2.0030303 

1. 9969696 

659.002 

661.002 

1.002 

It  is  evident  from  the  last  column  of  Table  VU.I 
that  the  lateral  freedom  decreases  linearly  with  the  a  - 
crease  in  time  lag  between  the  two  detectors.  Note  also 
that  the  lag  must  be  less  than  .25  seconds  for  the  lateral 


freedom  to  become  less  than  the  control  limits:  set  for 
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the  system. 


3.10  SAMPLER  AND  ZERO  ORDER  HOLD 


The  train  of  pulses  generated  by  the  detector 
activates  a  gating  circuit  to  the  input  buffer  register 
of  the  computer.  The  circuitry  is  represented  by  the 
block  diagram  of  Fig.  3*10. 


Output  of 
Detector 


rShaft  position  Encoder 

rTT  ?  Y 


plnput  Buffer  Register  |, 

FIGURE  3.10  BLOCK  DIAGRAM  OF  THE  INPUT  CIRCUITRY 

With  this  arrangement,  a  degree  of  frequency 
modulation  can  be  obtained  because  the  computer  program 
can  have  dummy  operations  to  make  the  time  for  successive 
samplings  equal „  The  computer  will  always  work  with 

the  latest  data  in  the  13R. 

Since  the  buffer  register  contains  information  that 

was  collected  some  time  prior  to  being  entered  into  the 
computer  for  processing,  and  since  the  processing  takes  a 
finite  time,  there  is  a  definite  time  delay  of  b  seconds 
before  the  signal  to  the  circuitry  constituting  the  guidance 
system  on  the  tractor,  is  transmitted.  Thus  the  computer 

imparts  a  time  delay  described  in  the  Laplace  domain  by 

-  -ST 
e 


«  0  * 


For  analysis  purposes,  the  time  lag  caused  by  the 

detector,  interface  equipment ,  and  the  computer  will  be 

bST 


dealt  with  as  a  lumped  transportation  lag  e 


in  -which 
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b  will  be  the  fraction  of  the  sampling  period  that  the 
time  lag  constitutes. 


The  pulses  generated  by  the  detector  as  its  drive 
system  oscillates  the  detector  field  across  the  path 
of  the  tractor,  are  to  be  used  as  carrier  pulses  to 
activate  gates  to  the  I3R  in  the  interface  hardware 
of  the  computer.  The  train  of  pulses  therefore,  ac¬ 
tivates  the  sampler  gates. 

To  simplify  the  analysis,  the  sampling  circuitry  may 
be  replaced  by  an  "ideal11  sampler  whose  output  signal  e*(t) 
contains  a  train  of  impulses.  This  approximation  can  be 
justified  if  the  pulse-width  of  the  sampler  output  is  very 
small  when  compared  with  the  dominant  time  constant  of  the 
continuous  signal  e(t),  and  to  the  sampling  period  T. 


e  ( t ) 


t) 


FIGURE  3*11  SAMPLER 


Assume  the  input  is  a  continuous  function  e(t)  as  in 
Fig.  3.11,  the  output  e*(t)  of  the  sampler  is  a  train  of 
impulses  and  is  related  to  e(t)  by 

e*(t)  =  e(t)6T(t)  3.10ol 


where 

6  T( t)  =  E {t -nT )  3.10.2 

n=-°° 

in  which  (t-nT)  represents  an  impulse  of  unit  area  oc¬ 
curring  at  time  t~nT .  Hence 
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3.10.3 


,  V  .  ’  00  -f  OO 

e*lt)  -  e ( t )  E«(t-nT)  =  £  e ( nT )6  ( t-nT ) 

n  =  -°°  n=-°° 

Assuming  that  e(t)  is  zero  for  T<0,  Eq.  3*10.3  can  be 
written  as 

°°  . 

e*  ( t )  =  E  e  ( n*T  )<5  ( t  -nT )  3.10.4 


Taking  the  Laplace  transform  of  e'*(t) 


directly  from 


Sq  3.10.4,  yields 

oo 

E*(S)  =L  (  e*  ( t ))  =LE{  e(nT  )6  ( t-nT)  }  3.10.5 

E*(S)  =  ?e(nT)e“nTs 

n=0  3.10.6 


If  the  Laplace  transform  of  e(t)  is  a  rational  function, 
E*(S)  can  be  written  in  a  closed  form. 

i.e.,  if  T  is  the  sampling  period  and  e(t)~u(t),  a  unit 


step  function, 

(S)  =  ? e ~nT s  =  i  +e-lTs  +e“2Ts  +e~3Ts  .,.3.10.7 


E 


1.7T 


n  =  0 


f  or 


3.10.8 


If  the  information  input  is  to  be  reconstructed,  it 


must  be  reconstructed  from  the  information  contained  in  the 


pulse  sequence.  Using  the  usual  method  for  generating 
this  approximation,  the  power  series  expansion  of  e(t)  in 
the  interval  between  sampling  instants,  nT  and  (n+l)T,  the 
following  equation  results: 

en(t)  =  e(nT)+e* (nT) ( t-nT)+e ' 1 (nT) ( t~nT)+. . . 

3.10.9 

where 

©n  ( t )  =  e(t)  for  nT<_t-<(n+l )T  3.10.10 


(3D 


' 


./hen  only  the  first  term  of  the  power  series  given  by 


no,  3.10.9  is  used,  the  extrapolatcr  is  referred  to  as 
a  zero-order  hold.  Then  Eq  3.10,9  can  be  written  as 

3.10.11 

The  impulse  response  of  a  zero-order  hold  circuit 
is  shown  in  Fig.  3*11 


en  ( t )  -  e  (  nT ) 


r 


r 


— .  t 


( a )  (  b ) 

FIGURE  3.12(a)  UNIT  IMPULSE  INPUT  TO  z-o-h;  (b)  II] 
PULSE  RESPONSE  OF  Z-o-h  TO  A  UNIT  STEP  INPUT 


The  transfer  function  of  the  hold  circuit  is 


Gho  ^ s '  “ 


-ST 


3. 10 o 12 


S 


The  combined  transfer  function  of  the  sampler  and  z-o-h 


elements  when  cascaded  is  obtained  from  the  following: 
S*(S)Gh0(S)  =_i —  rk£l  =  1 

1  ^  Q  Q 

1.  O  O  O 


3.10,13 


Thus  the  output  of  the  cascaded  elements  is  seen  to  be  a 
unit  step. 

A  mathematical  description  of  the  sampler  and  the 
z-o-h  have  been  developed  but  the  input  and  output  signals 
have  not  been  described.  The  origin  of  the  input  and 
the  destination  of  the  output  will  be  examined  in  the 
ensuing  sections 0 
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3.11  SHAFT  POSITION  ENCODER 

A  shaft  position  encoder  was  mentioned  earlier 
without  any  description  of  its  input-output  characteris¬ 
tics  other  than  that  its  output  is  a  binary  indication  of 
the  position  of  the  encoder  shaft.  The  Nordan  Ketay  shaft 
position  encoder  provides  in  binary  format,  a  digital  out¬ 
put  corresponding  to  the  position  of  the  encoder  shaft. 

The  position  of  the  detector  is  related  to  the  encoder  shaft 
position  through  a  gear  train  or  alternately  the  angular 
rotation  of  the  detector  can  be  related  to  the  number  of 
pulses  of  the  least  significant  bit.  The  former  method 
is  the  one  used  in  this  control  system. 

The  maximum  angular  velocity  of  the  shaft  position 
encoder  (hereafter  referred  to  as  the  SPE)  is  one  of  the 
factors  that  limit  the  input-output  ratio  of  the  gear  train 
between  the  detector  drive  system  and  the  SPE.  The  encoder 
produces  128  bits/turn  and  requires  64  turns  for  full 
count.  The  proposed  gear  train  is  not  restricted  to  any 
particular  ratio  because  the  conversion  of  the  output  to 
radians  is  accomplished,  internally  by  the  computer. 

Ratios  that  provide  a  low  angular  change  of  the  detector 
for  each  bit  of  the  SPE  are  desirable  provided  the  restric¬ 
tions  imposed  by  the  SPE  are  not  violated. 

The  only  restrictions  that  are  imposed  by  the  SPE 
are  (l)  the  SPE  not  be  operated  at  a  velocity  exceeding 

20Grpm.  and  (2)  practical  gear  ratios  be  used  (no  half 
teeth) 
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Fig*  3*13  is  a  sketch  of  the  driving  circuit  and  the 


first  three  stages  of  the  encoder  bit  positions. 


FIGURE  3*13  SPE  DRIVING  CIRCUIT  AND  THREE  STAGES 

TABLE  IX 


SPE  OUTPUT  RELATIONSHIP  TO  SHAFT  POSITION 


Input  position  0123  4  56789  10-11  12  13.... 


Bit  condition 


1 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

0 

•1 

A  •  ♦ 

2 

0 

0 

1 

1 

0 

0 

1 

1 

0 

0 

1 

1 

0 

0 »  . 

3 

0 

0 

0 

0 

1 

1 

1 

1 

0 

0 

0 

0 

1 

1 . . 

4 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

A 

5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.  , 

13  . * . .... 

The  bit  condition  table,  Table  IX,  confirms  that  the 
code  is  simple  binary;  thus  is  compatible  with  the  PDP-8 
computer  that  Is  referred  to  in  following  sections. 

Turning  at  the  maximum  velocity  of  200  rpm  the  encoder  pro 
duces  200x128/60  bits  per  second.  With  the  detector  scan¬ 
ning  9°/sec,  there  will  be  the  equivalent  of 

0°x60x60x. 000292/200x128  radians  per  bit. 

The  field  width  0  is  referred  to  again  in  later  sections. 
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3.12  GENERAL  DESCRIPTION  OF  THE  PDP-8  COMPUTER 

Analog  and  digital  elements  may  be  used  separately 


or  in  combination  with  one  another  in  automatic  control 


systemso  Likewise,  analog  and  digital  computers  can  be 
made  compatible.  An  investigation  by  the  author,  revealed 
the  following  advantages  in  using  digital  computers: 

(1)  as  much  precision  as  may  be  required  in  any  par¬ 
ticular  application  may  be  attained. 

(2)  improved  reliability  because  of  relative  freedom 
from  problems  encountered  in  electromechanical  or 
electronic  analog  computers,  namely  precision 

of  components,  wear,  and  stability. 

(3)  relative  insensitivity  to  mechanical  or  electri- 

'  cal  noise,  which  limits  the  complexity  of  problems 
solvable  on  an  analog  computer. 

(4)  facilities  exist  for  storage  of  large  quantities 
of  data  that  may  be  required  in  complex  problems-. 

( 5)  they  have  the  ability  to  perform  logical  as  well 
as  linear  or  nonlinear  mathematical  operations. 

(6)  they  have  flexibility  derived  from  their  logical 
capacity  and  the  fact  that  their  operations  are 
controlled  by  a  program  rather  than  specifically 
designed  or  interconnected  physical  components. 

The  flow  of  information  through  a  digital  computer 
system  is  illustrated  by  the  arrows  in  the  block  diagram 
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In  general  purpose  computers,  time  is  consumed  for 
ini  or  mat  ion  acquisit  ion ,  obtaining  access  to  data  and 
instructions  in  storage.  Also  in  the  arithmetic  operation 
considerable  amount  of  redundant  computing  must  be  done 
even  ix  the  values  of  the  input  information  change  very 
little  from  one  computation  to  the  next.  This  is  par¬ 
tially  oif set  oy  the  ability  of  the  computer  to  do  compu¬ 
tations  in  large  increments,  and  therefore,  the  sampling 

rate  in  a  control  sytem  can  be  varied  to  obtain  optimum 
stability. 


Precision  obtainable  in  a.  general  purpose  machine,  is 
governed  largely  by  the  numerical  method  and  program  selec¬ 
ted.  To  minimize  errors,  the  forms  of  the  equations  chosen 
must  be  such  that  they  adequately  represent  the  system’s 
characteristics ,  and  are  suited  to  the  characteristics  of 
the  computer. 


The  digital  computer,  for  which  the  following  work  is 
designed,  is  the  PDP-8  computer  manufactured  by  Digital 
Equipment  Corporation.  Minor  alterations  are  necessary 
for  the  system  to  be  compatible  with  other  general  purpose 
computers . 

3*13  INTERFACE  HARDWARE 

The  PDP-8  is  a  parallel-transfer  machine  that  distri¬ 
butes  and  collects  data  in  bytes  of  up  to  twelve  bits.  In 
the  Programmed  Data  Transfer  mode,  all  programmed  data 
transfers  take  place  through  the  accumulator.  The 
computer  program  controls  the  transfer  of  information.  In- 
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put  data  arrive  at  the  AC  as  pulses  received  at  the  inter¬ 
face  connectors  from  bussed  outputs  of  peripheral  devices. 
Gating  circuits  of  the  program  selected  device  produce 
these  pulses.  Command  pulses  generated  by  the  device  flow 
to  the  input/output  skip  facility  (IOS)  to  sample  the  con¬ 
dition  of  the  I/O  device  flags. 

3.13.2  Requirements  Imposed  on  Peripheral  Equipment 

The  peripheral  equipment  must  have  (l)  the  ability  to 
sample  the  select  code  and  when  selected,  to  be  capable  of 
producing  sequential  JOT  command  pulses  (device  selector); 
(2)  gating  circuits  at  the  output  of  the  transmitting  regis¬ 
ter  capable  of  sampling  the  information  in  the  output  regis¬ 
ter  and  of  supplying  a  pulse  to  the  computer  input  bus  when 
triggered  by  a  command  pulse  from  the  DS ;  and  (3)  a  busy/ 
done  flag  (flip-flop)  and  gating  circuit  which  can  pulse 
the  computer  input/output  skip  bus  upon  command  from  the 
DS  when  the  flag  is  set  in  the  binary  1  state  to  indicate 
that  the  device  is  ready  to  transfer  another  byte  of  in- 
formation . 

3.13.3  Device  Selector  (DS)  (W103) 

Bits  3  through  8  of  the  IOT  instruction  serve  as  a  de¬ 
vice  select  code.  Each  DS  is  assigned  a  select  code  and  is 
enabled  only  when  the  assigned  code  is  present  in  the  Memory 
Buffer  (MB).  When  enabled,  a  DS  regenerates  IOP  pulses  as 
IOT  command  pulses  and  transmits  these  pulses  to  skip,  in¬ 
put  gates  within  the  device  and  to  the  processor  to  clear 
the  Accumulator  (AG).  Fig.  3.1.5  is  a  symbolic  logic  repre- 
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sentation  of  the  selector. 


LOGIC  DIAGRAM  OP  THE  DEVICE  SELECTOR 
3.13*4  R123  Diode  Gate 

This  module  contains  six  dual  input  MAND  gates  for 
negative  levels  and  is  primarily  used  for  transferring 
data  into  or  out  of  the  PDP-8  accumulator.  Standard  DEC 
negative  levels  of  0.4  psec* negative  pulses  such  as  those 

from  the  W103  DS  can  be  used  as  input  signals.  Pig.  3,16 
is  the  logical  representation  of  this  module. 


NOTES 

’  STROPE  PULSE  APUT  TO  TERMINALS  F,  M,  AND  T  AmiCH  ARE 
CONNECTED  IN  COMMON  WHEN  USED  AS  A  BUS  GATE 

2  DATA  BiT  INPUTS  TO  TERMINALS  0,  E,  K.  L,  R.  AND  S 

3  TWO  MODULES  AR£  REQUIRED  TO  STROBE  A  12-BiT  *0*0 


Diode  Gale  R123  Logic  Circuit 


FIGURE  3*16  DIODE  GATE- 

The  strobe  pulse  is  input  to  terminals  F,  M ,  and  T.  Data 
bit  inputs  are  connected  to  terminals  D,S,K,L,R,  and  S. 
Two  modules  are  required  to  strobe  a  12-bit  word  into  the 
accumulator. 
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3*l4  COMPUTER  PROGRAM 


The  computer  program  must  instruct  the  computer  to  do 
the  following  operations: 


(1)  In  logic j  the  program  must  distinguish  between  the 
two  different  input s0 

(P)  Arithmetically  (a)  convert  the  input  to  operational 


language ; 


(3) 

with  the 
(  ^n~Rn~l{ 


(b) 

convert 

the  inputs  to  radians; 

(o) 

calculat 

e  the  cotangents; 

(d) 

sum  the 

cotangents . 

In  logic 

(a) 

compare 

the  quotient  of  AB-ct*. 

o  L 

previous  quotient  for  extraneous  results;  i.e»,  if 

is  greater  than  some  limit  of  validity,  wipe  out  ct- . 

1/ 

(b)  if  Jqn-Rn~lJ  is  less  than  the  con¬ 


stant  of  validity  then  go  on  with  the  next  comparison; 


(c)  compare  the  quotient  qn  with  the 
constant  (the  control  constant  r(t))  Rn-r(t)  =  e(t). 

This  error  signal  must  be  transmitted  to  the  steering 
amplifier  and  converted  to  analog  form  if  e(t)^K<^-e  (t). 

The  operation  of  selecting  an  input  device  was  refer¬ 
red  to  in  section  3.13.  The  DS  circuitry  is  similar  to 
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a  sequential  multiplexer  switch  except  that  the  computer 
program  calls  for  each  input  ab  its  convenience . 

The  conversion  to  operational  language  is  a  small 
operation.  The  computer  works  in  the  two’s  complement  of 
the  binary  input.  Thus  the  ones  are  converted-  to  zeros, 
the  zeros  to  ones  and  1  is  added  to  the  result. 

The  conversion  to  radians  is  accomplished  by  a  simple 
multiplication.  The  input  is  multiplied  by  the  conversion 
factor  D  referred  to  earlier  on. 

The  cotangents  of  the  angles  can  be  calculated  by 
solving  the  power  series, 


cot(x)=  1  _  X  _  x3  _  2x5  _  xg_  ....  3.14.1 

x  3  45  945  4725 

For  example,  if  cot ( 1.0)  were  required,  the  computer  would 
be  programmed  to  compute 

Got ( 1  o 0 )  1  _..1_  _  1  _  2  ,  1  3.14,2 

3  "45  945  4725 

Since  it  is  desired  to  know  cot(l.O)  to  within  *0.00001, 
the  series  would  be  terminated  after  the  term  that  had  a 
value  less  than  0., 00001.  Another  way  of  interpreting  this 
result  is  that  the  polynomial  approximation,  Eq.  3ol4„2 
is  sufficiently  accurate  for  our  purpose  when  x-1. 


Notation:  cot*x  indicates  an  approximation  to  the  cotangent. 
Another  method  that  was  considered  for  obtaining  the 
cotangents  of  the  inputs  was  by'  interpolat ion.  The  case 
has  been  considered  where  the  function  cot(x)  is  stored  in 
the  computer  as  a  subroutine.  There  are  obvious  limitations 
in  the  method  described.  The  routine  might  require  too 
much  time  or  too  much  memory  space.  In  these  circumstances 
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it  is  convenient  to  record  values 
memory  in  the  form  of  a  function  t 
corresponding  to  a  set  of  specific 


of  f(x)  in  the  computer 
able  which  gives  the  f(x^) 
values  of  x^.  It  appears 


^ u  fiifct  glance  t na o  this  nas  great  limitations,  since  only 
those  values  of  f(x)  which  correspond  to  a  relatively  few 
values  of  are  in  the  memory  of  the  computer.  However*  by 
means  of  interpolation,  values  of  f(x)  corresponding  to 
f s  not  listed  in  the  table  can  be  approximated. 

This  method  would  certainly  be  rapid,  however  the  use 
of  storage  units  is  not  particularly  economical,  and  the  ac¬ 
curacy  obtainable  would  be  limited  by  the  length  of  the 
word  used  in  the  computer.  Because  of  these  limitations, 
a  contraction  of  the  power  series  approximation  seemed  to 
offer  the  better  alternative  method  for  obtaining  the  co¬ 
tangents.  The  following  algorithm  is  a  contraction  of  the 


power  series. 

Single  Precision  Cotangent 
Algorithm: 


(1) 

If  x<  0 , 

use  cot(x)  = 

-cot ( x) ; 

Assume 

x^O . 

(2) 

Write  x 

=  tt  /4(  q )  +  r 

where  q  is  an 

integer  and 

0<  r  <tt/4 

.  Let  q_0  =  q 

( mod  4) . 

(3) 

If  c10  = 

0  or  2  (i.e., 

,  octant  1  and 

3) ,  define 

rQ  =  r.  If  q0  =  1  or  3  (i.e,,  octant  2  and  4), 
define  r0  =  V4  -  r. 

(4)  Define  the  case  number  s  as  follows;  for  cot(x) 
s=l  for  qo=0,  or  s=0  for  q0~l,  or  s=3  for  q0=2, 
or  s=2  for  q0=3» 

(5)  Compute  the  factor  P  as  follows: 
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if  rQ^  2 


3.14.3 


F  =  -4— ±— 13.-9.46ro  -  313.11 
r0  -  104.46  +  939,33/r~ 

If  rD<  2'  then  F  =  1 

The  maximum  error  of  this  formula  is  10“'. 


•"14 


(6)  Now  the  answer  is  r0/F  for  s=0;  F/r0  for  s=l; 

~F/r0  for  s-2 ;  and  ~rG/F  for  s=3. 

The  sum  of  the  cotangents  obtained  from  the  above  al~ 
gorithm  is  arrived  at  by  simple  addition.  The  quotient 
of  the  sum  divided  into  the  distance  between  the  detectors 
is  the  vital  byte  of  information;  the  steering  action  is  a 
function  of  the  difference  between  AB/ct^  and  the  constant 
h  which  represents  the  desired  perpendicular  distance  be¬ 
tween  the  base  line  AB  and  the  particular  path  rr *  that 
the  tractor  is  on.  The  sum  of  the  cotangents,  ct.^,  when 
divided  into  the  distance  AB  yields  a  value  to  compare  with 


the  desired  perpendicular  distance  h  in  Fig.  3.9. 

Ine  firs  c  check  that  must  be  made  on  this  byte  of  infor¬ 
mation,  is  regarding  its  validity.  If  |AB/ctt-h|>E  where 
e  is  a  constant  that  represents  the  bounds  of  validity,  the 
sum  would  be  discarded.  If  [AB/ctj.-h|<£  the  sum  is  valid. 

If  |AB/ct£~h|<£  and  |AB/ >  jy^  where  y  represents  the 
maximum  deviation  of  the  control  path,  the  sign  of  y  would 

be  transmitted  to  the  steering  control  system  on  the  trac¬ 
tor. 


If  ]AB/ctfc~h,l<  y ,  then  the  previous  control  signal  would 
be  changed  to  zero  or  left  at  zero  if  it  was  zero. 

ine  loutine  that  follows  immediately,  is  the  one  that 
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A 


w HI  ylsld  the  cotangent  of  angles  between  0  and  'n /4. 
minor  alteration  extends  the  range  of  angles  from  0  to  ” /A> 
to  0  to  t/2.  The  extended  program  is  referred  to  later. 
This  routine  was  prepared  in  order  to  establish  the  time 


required  uo  calculate  the  cotangents  to  the 
cision  using  the  above  algorithm. 


required  pre~ 


0 

0 

INPUT— 1 3 
OUTPUT- 1 4 

0007 

5600 

*1000 

TLS 

1000 

4407 

START , 

JKS  I  7 

1001 

0013 

INPUT 

1002 

6240 

FPUT  A1 

1003 

0013 

INPUT 

1004 

6243 

FPUT  A 2 

1003 

001.3 

INPUT 

1006 

6246 

FPUT  B2 

1007 

0013 

INPUT 

1010 

6251 

FPUT  B3 

1011 

0013 

INPUT 

1012 

6265 

FPUT  Cl 

1013 

0000 

F£XT 

101 4 

4407 

HERE , 

JMS  I  7 

1015 

0013 

INPUT 

1016 

6254 

FPUT  R0 

1017 

5254 

FG'ET  R0 

1020 

3254 

FMPY  R0 

1021 

62  57 

FPUT  R02 

1022 

5251 

FGET  B3 

1023 

4257 

PD IV  R02 

1024 

2246 

FSUB  B2 

1025 

1257 

FADD  R02 

1026 

6262 

FPUT  TEMP 

1027 

524o 

FGET  A1 

1030 

3257 

FMPY  R02 

1031 

2243 

FSUB  A2 

1032 

4262 

FDIV  TEMP 

1033 

1265 

FADD  Cl 

1034 

4254 

FDIV  R0 

1035 

0014 

OUTPUT 

1036 

0000 

FEXT 

1037 

5214 

JMP  HERS 
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. 


1040 

0000 

Al , 

0 

104l 

0000 

0 

1042 

0000 

0 

1043 

0000 

A2, 

0 

1044 

0000 

0 

1045 

0000 

0 

1046 

0000 

B2, 

0 

104? 

0000 

0 

1050 

0000 

0 

1051 

0000 

B3, 

0 

1052 

0000 

0 

1053 

0000 

0 

1054 

0000 

R0, 

0 

1055 

0000 

0 

1056 

0000 

0 

1057 

0000 

R02, 

0 

1060 

0000 

0 

1061 

0000 

0 

1062 

0000 

TEMP  f 

0 

1063 

0000 

0 

1064 

0000 

0 

1065 

0000  ■ 

✓*1  + 

e-1 » 

0 

1066 

0000 

0 

1067 

0000 

0 

Al 

1040 

A  2 

1043 

32 

1046 

B3 

1051 

Cl 

1065 

HERE 

101 4 

INPUT 

0013 

OUTPUT  0014 

RO 

1054 

R02 

1057 

START 

1000 

TEMP 

1062 

The  time 

required  f 

or 

cotangent  of 

the  angular 

i 

two  c 

0 tan gent 

s  are  requi 

re 

in  lo 

gic,  the 

total  time 

r 

the  t 

ime  dete 

rmined  above. 

The  next 

routine  is 

d 

logic 

required  in  the  oper 

exception  of 

filling  in 

-t-  ^ . 

orl 

4-  4-  V 
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ROUTINE  #2 


INPUT  =  13 
OUTPUT  =  l4 


*7 

5600 

*1000 

START , 

TLS 

JMS  I  ? 

INPUT 

FPUT  A1 

INPUT 

FPUT  A 2 

INPUT 

FPUT  B2. 

INPUT 

FPUT  33 

INPUT 

INPUT 

FPUT  Cl 

INPUT 

■ 

FPUT  SC 

scale  factor 

INPUT 

FPUT  P02 

pi  over  2 

INPUT 

FPUT  P04 

pi  over  4 

HERS, 

CLA 

clear  accumulator 

1ST 

strobe  input  buffer 

CIA 

2{s  complement  of  content 

JMS  I  7 

Call  floating  point  pack. 

INPUT 

FPUT  R 

store  R 

FGET  R 

FMPY  SC 

mult  by  scale  factor 

FSU3  P04 

FAC  =  R-Pi/4 

FEXT 

CLA 

TAD  HORD 

add  high  order  mantissa 

SPA 

skip  with  positive  ac 

JMP  LOW 

R0  =  R 

SMA 

JMP  COT 

cot  R  =  1 

JMS  I  7 

FGET  P02 

FSUB  R 

FPUT  RO 

Rq  =  Pi/2  -  R 

LOW , 

FEXT 

CIA 

DCA 

JMP  COMP  -1 
JMS  I  7 

JMS  I  7 

FPUT  RO 

Ho  =  R 

COMP, 

FGET  RO 

FMPY  RO 

Square  RQ 
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FPUT  R02 
FGET  B3 
PD IV  R02 
FSUB  32 
FADD  R02 
FPUT  TEMPI 
FGET  A1 
FMPY  R02 
FSUB  A 2 
FD IV  TEMPI 
FADD  Cl 
FPUT  TEMP 2 
FEXT 

TAD  CHECK 
SZA 

JMP  OUT 
JMS  I  ? 
FGET  RO 
FD IV  TEMP 2 
OUTPUT 
FEXT 
CLA 

DC A  CHECK 
JMP  HERE 

OUT,  JMS  I  7 

FGET  TEMP 2 
FDIV  RO 
OUTPUT 
FEXT 

JMP  HERE. 
COT,  CIA 

JMS  I  7 
INPUT 

FPUT  TEMP 3 

FGET  TEMP 3 

OUTPUT 

FEXT 

CLA 

DCA  CHECK 
JMP  HERE 


store  denominator 
c  ompu t  e  nu me  ra  tor 

compute  and  store  F 

check  whether  R  was  less 
than  Pi/ A 
no 

cot  R  =  Rq/F 


yes 

cot  R  =  F/Rq 


cot  R  =  1 


( A6) 


make  the  delay  time  an  integral  multiple  of  the  sampling 
period  of  the  detector  system. 

Mxplanat ions  to  the  right  of  the  instructions  are 
to  clarify  the  purpose  of  the  instruct ions .  Making  the 
estimate  of  the  time  required  to  complete  the  operations 
designated,  on ' the  basis  of  1.5  microseconds  for  logic 
operations  and  21.0  and  36.5  microseconds  for  multipli¬ 
cation  and  division  respectively,  the  time  reqired  is 
less  than  0.25  seconds.  Thus  if  a  sampling  frequency 
of  4  samples  per  second  is  adequate,  the  computer  is 
well  able  to  meet  the  specification. 
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CHAPTER  IV 


GENERAL  ANALYSIS  OF  THE  CONTROL  SYSTEM 

For  analysis  purposes,  using  approximations  to  the 
transfer  functions  of  the  equipment  is  the  rule  rather 
than  the  exception.  In  the  last  chapter,  approximations 
to  the  transfer  functions  of  the  equipment  were  established. 
Now  the  procedure  is  to  reestablish  some  order  to  the  con¬ 
trol  system  so  that  the  analysis  can  be  carried  out.  The 
block  diagram  of  the  control  system,  Fig,  4.1,  is  the  first 
step  in  that  direction. 

The  sampler  is  a  switch  that  is  activated  by  a  pulse 

that  originates  in  the  detector  amplifier  circuitry.  The 
bST 

time  lag  e  >  is  the  lag  imposed  not  only  by  the  computer 
but  also  the  transmitter;  the  transmission  of  infrared  rays; 
the  delay  of  few  micro-second  duration  for  which  the  data 
await  transfer  to  the  accumulator,  in  the  input  buffer 
register;  and  the  time  lag  of  the  response  of  the  guidance 
system  on  the  tractor.  Although  lumping  the  time  lags  into 
one  is  not  always  a  particularly  good  approximation,  it  will 
serve  in  this  case  where  the  undamped  natural  frequency  of 
the  system  is  relatively  low.  Further,  the  simplification 
achieved  in  the  analysis,  warrants  the  loss  in  accuracy  in 
the  approximation  of  the  response  of  the  system  to  a 


disturbance  o 


. 


4.2  GAIN  LIMIT  OF  A  SECOND-ORDER  CONTROL  SYSTEM 

Ihe  block  diagram,  Fig  4.2,  is  a  representation  of  the 
transfer  functions  of  the  elements  seen  in  Fig.  4.1. 

Fig.  4.3  is  a  simplification  of  the  block  diagram  in  Fig.  4.2. 
Ine  essential  difference  is  the  transfer  of  the  transpor¬ 
tation  delay  to  a  point  after  the  comparator  and  replacing 
the  synchronized  samplers  by  a  single  sampler  after  the 
comparator.  This  does  not  alter  the  sense  of  the  system 
since  R(S)  is  constant  over  a  time  period  Tp  T,  where  T 
is  the  time  period  of  the  sampler. 

R(S 


FIGURE  4.2  BLOCK  DIAGRAM  OF  THE  ELEMENTS  PRESENT 

IN  FIGURE  4.1 


FIGURE  4.3  SIMPLIFIED  BLOCK  DIAGRAM  OF  FIGURE  4.2 
TJn  in  the  above  figures,  is  a  nonlinear  coefficient 
that  describes  the  output  of  the  relay  with  dead  zone, 
which  is  represented  in  Fig.  4.2  and  4.3  by  N.  The  linear 
forward  function  G(S)  is  the  approximation  of  the  transfer 
function  of  the  guidance  mechanism  on  the  tractor.  For 
this  general  analysis,  the  transfer  function  of  the  steering 
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mechanism  is  rj'mv2/s  (s+l ) .  The  unity  damping  is  provided 

by  the  gyrocompass  controlling  a  potentiometer  in  a  servo 
circuit. 


It  is  noted  that  the  transfer  function  of  the  detector 

crlve  system  does  not  appear  in  the  transfer  function  of 

the  above  system.  It  is  assumed,  for  analysis  purposes, 

that  the  two  samplers  are  synchronized  however. 

^.3  GAIN  LIMIT  OP  A  SAMPLED  DATA  CONTROL  SYSTEM  WITH 
ZERO-ORDER  HOLD 


Consider  the  second-order  sampled  data  control  sys¬ 
tem  in  Fig,  4.5, 


FIGURE  4.5  SAMPLED  DATA  SYSTEM  WITH  ZERO-ORDER  HOLD 
The  transfer  function  of  the  forward  path  is  GhoG(S). 
That  means 


C*(S)  =  GhoG*(S)E*(S) 

The  z-transform  of  C(s)  is  then 

C(Z)  =  AMPM  H(Z) 

1  +  HGh0G(Z) 

This  can  be  written  as 


4.3. 1 

4.3.2 


C  ( Z )  GV)  rvG  ( Z ) 

R(zT  s  I"+  hgj^gTzT 

which  is  recognized  to  be  the 
in  Fig.  4.3.  Adding  the  z-o~h 
a  sampler  gives  an  additional 


4.3.3 


z~ transform  of  the  system 
»  Gho»  t 0  the  circuit  with 
term  1 to  the  forward 

. s . 


transfer  function. 


(5D 


N  ow 


r*  (  o  \ 

U'no'j|.oj 


-St 


4.3.5 


s  +  1) 


S  ince 


C(Z) 


K(Z) 


'no  - 


r-(Z) 


4.3,6 


C(Z)  (1 


Z"1)  ,  Gh0G(S) 

’ '  r  1  +  HG  G 


R(Z) 
in  which 

Z-ho^(S)  ~  2HGho^(S)  =  ^ 
Ghn-(Z)  = 


(3) 


4.3.? 


K 


(1 


s^(s  -5-  1) 


,  - .  r  tz 

( l-z~l)Kf - — 

1 . L  ( z  . 

•  1)' 

(1 


z”1! 


"T)z  ) 


4.3*3 

'+.3.9 


^~G(Z) 


'-hov 


(  Z  —  1  )  (  Z  —  G  *' 

K(  (T~l*e~T  )  z  -Te“J+l-e~’^ )  4^  3  ^  q 

( z-l) ( z-e~T) 

The  characteristic  equation  obtained  by  setting  the  denomin¬ 
ator  of  Eq.  4. 3*6  equal  to  zero  is 


(z-l) (z-e~T)  +  K((T-l+e“T)z  -Te  +l-e“T)  =  0  4.3* 


11 


To  obtain  the  relationship  between  K  and  T,  they  are 
left  unspecified.  It  is  also  necessary  to  change  the  equation 

to  the  r-transform  in  order  to  apply  the  Routhian  technique. 
This  is  done  by  substituting  ( r+1 )/( r-1)  for  z  in  Eq.  4.3*11* 

2  _  (r-H)„e~T  K(T-l*e"T) (r+l) 


(r-1)  (r-1) 


( r-1 ) 
2 


-T 

•Te  -M 


0  4.3*12 


Multiplying  Eq.  4,3*12  by  (r-1)  yields 

2  r-r  2  -2  re  ~T+  2  e  “T+ K  ( T  - 1 + e  -T )  r  ~K  ( T  - 1+  e  -T )  -  (Te  -T  - 1+ e -T ) r 
+2K(Te”T-i+e-T)r  -K(Te-T_i+e-T)  =  0  4,3*13 
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Thi 


Rout'nian  array  of  Sqfc 
KT  ( 1 ~  e “^ ) 


4.3.13  is 

2 ( l+e“T) -K(Te~T~2+2e~T+T) 


r1  2~2e“?-5-2K(Te**T-l+e“T) 
r°  2  ( 1+  s  ~T  )  -X  ( T  e  ~T  -  2+  2  e  "T+T ) 


S  ince 


uiiere  are  no  negative  signs 


in 


the 


P  • 
x  . 


column 


of  the  array,  it  is  concluded,  that  there  are  no  negative 
roots  for  a  certain  range  of  values  of  K  and  Y.  The 
Table  X  shows  numerically  the  relationship  between  speci¬ 
fied  values  of  T  and  K  in  the  r-  and  r^  coefficients  of  the 


array.  The  data  is  displayed  graphically  in  Pig.  4*6. 
TABLE  X  MAXIMUM  GAIN  VALUE  FOR  SPECIFIED  T 
T  1  .9  .8  o7  .6  .5  .4  .3  *  2 


Kr1  2.39  2.67 

2.88 

3.25  3.70  4.37  5.88 

7.40  11.3 

Kr2  27  35 

52 

75 

FIGURE  4.6  GAIN  VERSUS  SAMPLING  PERIOD 


Now  suppose  the  hold  function  ( l-e”^3) »  is  delayed 

s 

by  for  example.  What  is  the  net  effect  on  the  gain? 


Gh0G(S)  -  (l-e-3T)e0-5sT(  k) 


s 


sTs+1) 


4.3.14 
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« 


Tne  characteristic  equation  1+HG^qG(Z)  —  0,  becomes 
(  z-1 )  (  z-e"1  )  +Kc^‘^jl(  (T”lve”'^)z  Te^+l-e”^)  -0 


4.3 


-  r 
•  JLO 


Conversion  to  the  r-transform  results  in  the  equation 

go.  5inKT(  l-e~T)r  +  r(2(  l+e“T )  +  e° •  5%(Te~T~2+2e“T+T )  ) 


+  2( l+e“T ) “e0  *  ^TK(Te“T-2+2e“T+T )»0 

4.3,17 

The  effect  of  the  transportation  delay,  e~°*^-\  in  the 
system  is  to  multiply  K  by  a  factor  of  1/1.6487  (without 
sacrificing  stability)  when  T  is  one  second.  Table  XI 
shows  numerically,  the  relationship  between  T  and  the  am¬ 
plification  factor  when  the  transportaion  delay  is  \  sec. 


The  data  is  graphically,  displayed  in  Fig.  4.5. 


TABLE  XI  MAXIMUM 
T  1 

F(eisT)  i e 28 
K-max  i  « 8  7 


K  FOR  SYSTEM  WITH 
.75  .50 

1.35  1.65 

2.22  2.65 


TRANSPORTATION  LAG 
.25 
2.71 
3 . 4o 


•  Z$  '$o  .*£  / 


FIGURE  4.7  GAIN  VERSUS  SAMPLING  PERIOD  IN 
A  SYSTEM  WITH  TRANSPORTATION  LAG. 

From  the  examination  of  the  system  with  transporta¬ 
tion  lag,  it  is  found  that  the  system  is  stable  for  certain 
limits  of  gain.  Since  the  gain  is  TmvZ/s{ s+1 ) ,  in  which  T 
is  the  steering  constant  and  v  is  the  velocity  of  the  trac- 


i 


tor,  seen  that  the  limits  of  gain  can  be  realized, 

oy  decreasing  either  the  steering  constant  or  the  veloci¬ 
ty  of  the  tractor.  Hence,  the  system  theoretically  works 
satisfactorily  within  limits  of  speed  and  turning  radius 
that  are  practical  in  field  operations. 
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CHAPTER  V 

RESPONSE  DETERMINATION 

It  was  established  that  the  control  system  may  be 
designed  to  be  stable  for  some  operating  conditions.  The 
operating  conditions  that  the  system  would  be  subjected  to 
in  practice,  however,  may  be  nowhere  near  the  stable  range. 
Operating  in  the  realm  of  possibility  that  the  coefficients 
that  are  derived  theoretically,  are  also  correct  for  a  prac¬ 
tical  system,  observations  of  the  response  of  the  control 
system  that  has  only  the  initial  conditions  as  disturbances, 
are  of  interest.  It  is  recognized  that  in  practice  there 
is  virtually  no  linear  system;  however,  approximations  of 
practical  operating  conditions  by  linear  functions,  will 
lead  to  reasonable  conclusions  regarding  the  performance  of 
the  system  provided  the  limits  of  operation  do  not  deviate 
from  the  region  of  approximation. 

With  the  above  realization,  the  operating  conditions 
listed  below  are  derived. 

Operating  Conditions: 

(1)  i  mile  square  plot,  Fig.  5*1 

(2)  operating  at  an  extreme  point 

(3)  it  is  desirable  to  leave  a  headland  of  20  ft. 

(4)  there  is  unobstructed  view  to  all  parts  of  the  plot 

(5)  the  tractor  is  located  at  point  P,  20  ft.  in 
from  an  edge  of  the  plot 


(6)  there  are  observation  points  at  A  and  B,  the 
end  points  of  the  reference  line,  1320*  apart 
and  1320*  from  the  desired  path  rr* 

(?)  the  operating  limits  are  1319  ’< h<132l 1 

(8)  the  initial  conditions  are  x^-1321  and  X2=0 

(9)  the  modulated  sampling  period  T,  is  1  second 

(10)  the  turning  radius  of  the  tractor  is  20  ft. 

(11)  the  velocity  of  the  tractor  is  8  ft. /sec 

(12)  the  gyrocompass  feedback  is  -1 

(13)  from  (10)  and  (ll)  the  steering  constant  Tm 
is  found  to  be  c05  radians/foot. 

(14)  the  lumped  time  lag  is  .25  sec. 


The  transfer  function  of  the 
tractor  is  G^(S)  -  64/s(s+l).  The 
the  relay  with  dead  zone  are  +1.0. 

ft -a 

The  control  system  is  represented 


guidance  system  on  the 
switching  limits  of 
R(tQ)  is  a  unit  step, 
in  the  block  diagram  of 


Fig.  5.2 


FIGURE  5.2  BLOCK  DIAGRAM  OF  THE  CONTROL  SYSTEM 
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5.1.2  AM A LYSIS  OF  THE  SAMPLED -DATA  SYSTEM  WITH  NON LINEARITY 
BY  THE  STATE  TRANSITION  TECHNIQUE  WITH  VARIABLE  GAIN  CONCEPT 

The  variable . gain  concept  with  the  state  transition 
method,  permits  a  systematic  synthesis  procedure  for  dis¬ 
crete-data  control  systems.  The  variable  gain  concept  is 
a  method  of  representing  a  nonlinearity  by  a  variable  gain 
amplifier  with  gain  K(kT)  =  during  the  sampling  period 
kT  < t £  (k+l)T.  Use  is  made  of  the  transition  flow  graph 
and  equations,  thus  the  method  is  essentially  numeric. 

Kv  is  defined  as  Kv  =  TTpf£'.tJ.  5*1*1 

K  h(kT+) 

where  m  is  the  output  of  the  nonlinear  element  and  h  is 
the  input  to  the  nonlinear  element  from  the  z-o-h."  An 
element  with  similar  characteristics  to  the  one  used  by 
MacHardy-^*  in  the  continuous  data  system  will  be  con¬ 
sidered.  The  nonlinearity  is  a  relay 

type  function  with  dead  zone  and  with  its  input-output 
characteristics  described  by 


m(t)  = 


From  the  transition  flow  graph,  Fig.  5*3 »  of  the 
nonlinear  sampled- data  control  system  in  Fig.  5*2,  the 
transition  equations  are  obtained  in  the  following  vector 
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FIGURE  5.3  SIGNAL  FLOW  GRAPH  OF  THE  SAMPLED -DATA 
CONTROL  SYSTEM 
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matrix  form: 

x(t)-  K((t-tQ)  ,K(t0))xi(t0)+  N((t-tQ)  ,K(fc*))r(t0) 

5*1*3 

when  t"-  ( k+ 1 )  T  and  tQ=kT  ,  Eq ,  5*1.3  becomes 

x ( (  k+ 1 )  T )  -  M(T,Kk)Xi(kT)  +  N(T,Kk)r(kT)  5.1.4 

The  values  of  are  given  by  Eq.  5.1.1  in  which 
h(kT+)  is  determined  from 

h(kT+)  -  r(kT)  -  xx(kT)  5.1.5 

and  m(kT*t)  is  determined  from  the  nonlinear  characteris¬ 
tics  once  h(kT+)  is  known.  After  the  values  of  Kv 
during  each  sampling  period  are  determined,  Eq.  5.1,4  gives 
the  solution  of  the  state  variables  at  the  sampling 
instants . 


For  the  state  transition  flow  graph  of  the  system 
drawn  in  Fig,  5. 3s  the  transition  equations  are 

xx(  (k+l)T)=(  1-Kk *v2 (T~l+e~™ )  )xx(kT)  +v(  l-e~~  ) x2( kT ) 

+Kkv2 ( T-l+e“T ) r( kT )  5.1.6 

x2 ( ( kU ) T ) = -Kvv ( i -e"T ) xx ( kT ) + e “Tx2 ( kT ) 

+Kkv  ( l-e'*1 )  r( kT )  5.1.7 

If  the  initial  conditions,  T=1  sec.,  v=8ft/sec,  and  r(kT)=0 

are  substituted  in  Eq.  5*1.6  and  5*1.?»  they  become 

x i  ( ( k+ 1 ) T ) ~ ( 1—64*0.3 68 Kv ) x i ( kT )  +  8*0.632x2(kT) 

1  ‘w  5.1.8 

x2.(  (k+l)T)=(-8*0.632Kk)xx(kT)  +  0.368x2(kT)  5.1.9 

In  the  analysis  of  the  control  system  the  time  delay 
must  also  be  considered.  Recall  that  the  time  delay  is  the 
period  between  the  time  the  position  of  the  tractor  is  sam¬ 
pled  and  the  time  the  signal  reaches  the  guidance  system 
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on  the  tractor.  It  was  determined  that  the  lag  can  be 
limited  to  approximately  0.25  seconds  in  Chapter  IV. 
io  cope  with  the  delay  of  0.25  seconds  in  the  analysis, 
the  was  calculated  on  the  basis  of  the  conditions 
of  X1  at  the  instant  when  t=3/4(T).  The  new  Kk  was 
then  substituted  in  the  equations  5.1.8  and  5.1. 9, 
and  the  new  X^  and  X2  were  calculated  on  the  basis  of 
the  values  of  Xj  and  X2  dtermined  at  the  end  of  the 
last  sampling  period  T. 

The  set  of  equations,  Eq .  5.1.8  and  5*1.9.  were 
programmed  on  the  APL/360  and  solutions  for  the  first 
20  intervals  or  sampling  periods  were  calculated.  The 
routines  used  on  the  APL/360  computer  and  the  results 
produced  are  found  in  the-  Appendices.  Fig.  5. 4,5. 5,  and 
5.6  are  plots  of  the  displacement  as  a  function  of  time. 
5.2  DISCUSSION 

Since  there  is  no  reason  to  anticipate  occurrence  of 
frequencies  less  than  the  sampling  frequency,  the  values 
of  Xi  can  be  connected  by  smooth  curves  as  was  done  in 
Figures  5*^»  5*5.  and  5*6.  The  smooth  curves  produced, 
are  a  reasonable  estimation  of  the  performance  of  the 
system  between  the  sampled  instants.  The  equations 
were  not  solved  for  more  than  20  sampling  periods  because 
information  that  could  be  obtained  by  carrying  the  solu¬ 
tions  out  to  infinity  would  be  questionable  on  the 
basis  of  the  approximations  that  were  made  to  obtain  the 
mathematical  description  of  the  system. 
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Solutions  to  the  equations  were  obtained  for  only 
one  set  of  initial  conditions;  a  lateral  displacement, 

,  of  one  foot  and  zero  lateral  velocity,  X2*  Initial 
conditions  other  than  those  above  would  lead  to  different 
transient  solutions  of  to  different  periodic  solutions, 
according  to  Kuol^. 

CONCLUSIONS 

Based  on  observations  of  the  solutions  to  the  set  of 
Eq’s,  5*1c8  and  5« 1*9,  it  is  concluded  that: 

(1)  Operating  at  8  ft/sec  and  a  sampling  frequency  of  once 
per  second,  the  control  system  has  an  unstable  transient 
response.  However  by  increasing  the  sampling  frequency  to 
2.5  times  per  second,  the  transient  solution  can  be  made 
stable.  Thus  increasing  the  sampling  frequency  leads  to 

a  more  stable  transient  response. 

(2)  Besides  leading  to  a  more  stable  transient  solution, 
an  increased  sampling  frequency  results  in  lower  overshoot, 
smaller  steady  state  error,  and  shorter  settling  time. 

(3)  Wider  dead,  zone  limits  lead  to  larger  steady  state  error 
but  a  more  stable  transient  response.  For  example,  the  system 
when  operated  with  a  dead  zone  limit  of  ±2  feet,  is  stable 

at  a  velocity  of  8  feet/second,  and  sampling  period  1  sec. 

(4)  Longer  turning  radii  lead  to  more  stable  response,  lower 
natural  frequencies,  lower  overshoots,  and  to  larger  steady 
state  errors. 

(5)  By  increasing  the  demands  on  the  detection  system  so  that 

the  position  of  the  tractor  can  be  sampled  every  0.4  seconds, 

(64) 


- 

*  . 


the  tractor  has  a  stable  transient  response.  However,  to 
increase  che  sampling  frequency,  it  is  necessary  to  reduce 
the  field  scanned  by  the  detector  in  order  to  maintain  the 
same  relative  position  of  the  tractor  exhaust  muffler  in  the 
detector  field  when  the  $0%  level  of  response  is  reached. 

For  example :  a  field  width  of  12-i°  can  be  scanned  when  the 
sampling  frequency  is  once  per  second,  the  region  that  can 
be  scanned  is  reduced  in  width  to  5°  (due  to  the  maximum 
velocity  of  scan  of  ,22  radians  per  second). 

ADDENDUM 

The  control  of  a  tractor  by  the  method  described  in 
this  thesis  is  possible;  how  practical,  however,  remains 
to  be  judged  in  each  application  under  consideration. 

Field  operations  with  wide  implements  that  require  control 
in  overlap  of  less  than  one  foot,  can  be  carried  out  by  the 
automatically  controlled  tractor*  but  between  the  rows  culti¬ 
vation  (for  example)  requires  greater  precision  than  that 
designed  into  the  system  described  in  this  thesis. 
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